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Abstract 

Can cloud screening of an aerosol data set, affect the 
aerosol optical thickness (AOT) climatology? Aerosols, 
humidity and clouds are correlated. Therefore, rigorous 
cloud screening can systematically bias towards less cloudy 
conditions, underestimating the average AOT. Here, using 
AERONET data we show that systematic rejection of 
variable atmospheric optical conditions can generate such 
bias in the average AOT. Therefore we recommend (1) to 
introduce more powerful spectral variability cloud 
screening and (2) to change the philosophy behind present 
aerosol climatologies: Instead of systematically rejecting 
all cloud contaminations, we suggest to intentionally allow 
the presence of cloud contamination, estimate the statistical 
impact of the contamination and correct for it. The 
analysis, applied to 10 AERONET stations with ~ 4 years 
of data, shows almost no change for Rome (Italy), but up to 
a change in AOT of 0.12 in Beijing (PRC). Similar 
technique may be explored for satellite analysis, e.g. 
MODIS. 

Introduction 

Evaluation of chemical transport models and initiation of 
climate models is based on average aerosol properties 
measured from ground based or satellite platforms, over a 
given time and space interval (Chin et al., 2002, Menon et 
al., 2002). AERONET (http://aeronet.gsfc.nasa.gov/) cloud 
screening is based on variability of the 1 minute interval 
triplet measurements and -15 minute interval 
measurements (Smirnov et al., 2000). The basic 
assumption is that clouds vary more than aerosols and a 
given threshold of variability can separate clouds from 
aerosol. A similar technique is used for MODIS, though the 
MODIS algorithm, by using two-dimensional spatial 
variability method is more powerful (Martins et al., 2002). 
Kaufman et al. (2005a) used a new spectral variability 
cloud-screening algorithm (SVA) of AERONET optical 
thickness data (Holben et al., 1998; 2001). Application to 
one month of data collected in Lille, France indicates that 
the present LI .5 AERONET cloud screening may not retain 
variable aerosol and reject it as clouds. Aerosols can vary 
due to variability in humidity, the presence of near by 
sources or cloud processing. Due to the correlation between 
the aerosols and cloud cover (Chou et al., 2002; Sekiguchi 
et al 2003; Kaufman et al., 2005b; Koren et al., 2005), time 
averaged values may be influenced by the particular use of 
cloud screening. 

Spectral variability cloud screening algorithm 

The SVA (Kaufman et al., 2005a) can represent better the 
aerosol climatology in case of aerosol with spectrally 
varying optical thickness. The physical principle behind the 
SVA is shown in Fig. 1 . Quick observation shows that the 
main difference between the clouds and the aerosol in the 
MODIS image is the difference in color. We can also see 
that while the clouds are more variable than the aerosol, the 


heavy pollution aerosol is also highly variable. Elimination 
of the variable aerosol means elimination of the most 
concentrated aerosol. The SVA is designed to screen as 
clouds only measurements with significant spectrally 
neutral variability. Plot of the AERONET data using the 
SVA and L1.5 cloud screenings (Fig. 1) shows under 
representation of the hazy conditions in the LI. 5 data. 
There is no indication in the time dependence of the 
Angstrom exponent to indicate that some of the points are 
cloud contaminated. 
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Fig. 1: Top- MODIS image (July 20, 2005) of pollution in 
Beijing region (blue circle) observed from the TERRA satellite. 
The image shows the highly variable white clouds and still 
variable blue pollution aerosol with aerosol optical thickness 
as high as 5.0. The MODIS image, downloaded from the 
AERONET web site are enhanced for better spectral contrast. 
Bottom- AERONET AOTs using SVA (diamonds) and LI. 5 
(circles) for this day. And Angstrom exponent (triangles). 


Another demonstration of the need of spectral screening in 
more controlled conditions is shown in Fig. 2. The 
measurements represent thin smoke plumes generated by 



numerous small fires in the Lag-Ba-Omer holiday in Israel. 
The spectrally sensitive SVA detects most of the smoke 
plumes while AERONET LI. 5 rejected 2/3 of them as 
clouds. We can see in the figure that the rejected data, 
though variable, have much stronger variability in the blue 
spectral range than in the NIR, resulting in similar 
Angstrom exponent and indicating variable aerosol rather 
than clouds. Similar rejection of smoke plumes was found 
by O’Neill, etal. [2003]. 
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Fig. 2: Time variability of the spectral optical thickness 
and Angstrom exponent (triangles) in Nes Ziona, Israel 
during the holiday of Lag-Baomer when fires are set and 
result in numerous thin smoke plumes. SVA (filled squares) 
detected most of the measurements as cloud free while less 
than 1/3 were detected as cloud free by AERONET LI. 5. 
The LI. 5 average AOT is about 25% smaller. 

The SVA used in Figs 1-2 is applied to the level 1.0 
AERONET data (that already have some preliminary high 
frequency variability screening) using the following criteria 
(after Kaufman et al., 2005a): 

Triplet screening (variation over seconds) using spectral 
neutral variability: 

It is a cloud if: 6t cloud > 0.005+0.02x 675 (1) 

With &t cloud = 5x 870 - bt 44 V 71 7x 440 ) 

where inf is the spectral triplet value for wavelength X, x x 
is the measured optical thickness and 6x dou<i is the estimated 
variability of the cloud optical thickness. The threshold 
dependence on the AOT in eq. 1 represents the effect of 
humidity on the spectral dependence of the AOT. The 
cloud variability 5x doud (eq. 1) is the spectrally neutral 
component of the triplet variability. Lets take an example 
of aerosol with A=l, and with the triplet variability 
originating only from thin cloud that is present only part of 
the time. Then •t 87 °~0.5T ,MO and 6x 87Q -bT i4 °. Calculating 
&T cl °“ d we get half of the actual cloud optical thickness in 
this case. In another example, of pure aerosol with A=1 but 
variable in time, we shall get x 87 ”~0.5x 440 and 6x 87 °~0.56x 440 
resulting in 8x d ° ud =0. 


Adjacent screening (variation over minutes): 

It is a cloud if: Ax cloud > 0.007 5+0.03x 675 (2) 

with Ax cloud =Ax E70 -Ax 440 (x 870 /x 440 ) 

where Ax’ 1 is the maximum difference between the current 
AOT and the next or previous one. Ax do '“ 1 is the estimated 
variability of the cloud optical thickness. Only data with 
Angstrom exponent, A(440-870) > 0.3 are analyzed here 
(A<0.3 represents clouds or pure dust conditions). If the 
screened value is found later to be surrounded by values 
declared as cloud contaminated then the value is 
eliminated. To be consistent with LI. 5 data we screen out 
measurements for solar zenith angle > 78.5°. 

Application of the SVA algorithm is shown for 4 locations 
in Fig. 3 and compared with the AERONET LI. 5 
algorithm. The SVA generates larger AOTs in Beijing and 
ISPRA for A<1. These were found to be the heavy 
pollution conditions indicated in Fig. 1. The density of 
measurements increases all along the A axis, and is most 
pronounced in Beijing for A~1 .3 and in Alta Floresta in the 
presence of biomass burning smoke for A-2.0. 

Aerosol climatology 

The density of measurements, the average AOT, and A 
depend on the threshold of the SVA (Fig. 3). Thus even 
though SVA allows variable aerosol to be considered as 
aerosol rather than clouds, in the presence of cloud 
contamination the variable aerosol is still screened out with 
the clouds. Is there a way to avoid this trap of heavy 
aerosol being thrown away with the clouds? In Fig. 4 and 
Table 1, we explore variations of the average AOT and the 
average A with variations of the cloud screening 
thresholds. The results are plotted as a junction of the 
estimate of the cloud contamination. To estimate the cloud 
contamination we calculate the average cloud variability 
for all the measurements with positive 5x doud or Ax doud . The 
cloud variability is then converted into the cloud optical 
thickness using the relationship shown in Fig 4 bottom 
panel. Here the lidar measured variability of the cloud 
optical thickness is plotted as a function of the actual cloud 
optical thickness. This relationship basically suggests that 
the cloud optical thickness is on average twice the 
variability of the cloud optical thickness across 1 minute. 

As expected, relaxation in the SVA threshold increases the 
average AOT and decreases A, a clear evidence of the 
cloud contamination. However, correcting for the cloud 
contamination corrects for the change in A, leaving it 
practically constant and still leaves in some cases a higher 
value of the AOT. Note that the minimum value of the 
AOT indicates an optimum selection of the cloud screening 
threshold with minimum contamination. The results of the 
process for 10 AERONET sites are summarized in Table 1. 
For some of the locations (Rome, Kanpur, Mongu) the 
SVA did not change significantly the average AOT but 
increased the data rate by 4-20%. For other sites SVA 
increased the average AOT by 0.02 or 9% for Alta Floresta 
up to 0.08 or 20% for Beijing. Relaxing the SVA threshold 
and correcting for the cloud contamination still did not 
make a significant difference for Rome, Kanpur & Mongu 
but increased the average AOT for other site with a 
maximum increase of 0.12 for Beijing. The large difference 
in Beijing between the LI. 5 and the SAV algorithm can be 
traced back to Fig. 1. The heavy aged pollution in Beijing 
with AOT values as high as 3.0 are sub-sampled in the LI 5 
algorithm better represented by the SAV algorithm. 
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Fig. 3: Aerosol optical thickness (light color solid lines) and density of measurements (strong color dashed lines), dn/dA, as a 
function of the Angstrom exponent, A. The data, from 2001-2005, are sorted by A and averaged in groups of 100. ), dn/dA is 
the number of measurements, n, per unit A. Black/gray - AERONET level 1 .5 data; red/orange - spectral variability 
algorithm (SVA); green - SVA for thresholds reduced by 33%; blue/aqua - SVA with thresholds increased by 50%. The high 
AOTs with A~1 .0 in Beijing, China and Ispra, Italy (Melin and Zibordi, 2005) are due to accumulation of pollution over 
China and under the Alps respectively. Low A values (-0.5) are due to desert dust. Alta Floresta, Brazil represent an area with 
concentrated biomass burning aerosol. Rome represents mostly local pollution with dust intrusions for lower values of A. 

Table 1: The average aerosol optical thickness (AOT), Angstrom exponent (A), and percent of measurements after cloud 
screening (%) with A>0.3. Four estimate are made: (1) Average after SVA screening, (2) SVA corrected for residual cloud 
contamination, (3) no cloud screening except of A>0.3, and (4) the AERONET LI .5 data. Comparison among the values of A 
among the methods for a given site is used as a sanity check regarding residual cloud contamination. 
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Discussion 

The SVA method and the inclusion of all the data 
practically with no further cloud screening of the 
AERONET level 1.0 data (for A>0.3) in the analysis of 
long-term averages are based on existence of spectral 
difference between aerosol and clouds. Therefore we 
limited ourselves to A>0.3 and were not able to include 
pure dust such as observed in Capo Verde with A<0.3 in 
our analysis. In Fig. 5 we analyze the statistics of 
occurrence of A>0.3. Except of Capo Verde, the restriction 
of A>0.3 retains 88-98% of the data. In Kanpur the 
restriction of A>0.3 will exclude most of the dust in the 
spring time. Correction of the cloud contamination is based 
on a factor of 2 found to relate the cloud optical thickness 
and the variability of the cloud optical thickness across the 


1 minute of AERONET observations. This was measured 
by lidar on stratocumuli in 3 different locations (maritime, 
rural, and semi-urban) and altitudes, 600, 2500 and 1000 m, 
respectively. This relationship should be established in 
more numerous cloud conditions. However, the fact that the 
corrected average value of A are practically independent of 
the threshold of cloud screening indicates that the 
correction is performing well. 

Summary 

The correlation among aerosol, the humidity field and 
clouds introduces errors into aerosol climatology derived 
using rigorous stringent cloud screening. Using AERONET 
data we show that screening algorithm that uses spectral 
variability is able to discriminate better between clouds and 






















































aerosol, but probably the best climatology can be achieved 
using a relaxed cloud screening and estimating statistically 
the cloud contamination, subtracting it later from the 
climatology. We showed that using this technique we 
derive higher aerosol average AOT in some locations while 
keeping the same values in others. Much larger fraction of 
the data are used in the process, while still maintaining 
fairly constant Angstrom exponent. Similar techniques may 
be also applied to satellite data of aerosol, e.g. MODIS 
(Renter et al 2005). 
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Fig. 4: Top - Aerosol optical thickness (circles) and 
Angstrom exponent (diamonds) calculated for a range of 
SVA cloud screening thresholds and plotted as a function of 
the average residual cloud contamination. Full symbols — 
corrected for clouds. Open symbols - no correction for 
clouds. Bottom - Lidar measurements of the 1 minute 
variability in the cloud optical thickness (ACOT) as a 
function of the actual COT (for ACOT<l). Different 
symbols are for 4 different days and locations of 
measurements. The power law fit is for all the data. 
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Fig. 5: Bottom: Cumulative histograms of the data as a 

function of the Angstrom exponent for several locations. 

The cut off of A=0.3 except of Capo Verde retains 88-98°% 

of the data. 
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